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Abstract 

The connection between discrete and continuous dynamical systems 
through the unit-time map has shown a significant role in bifurcation the- 
ory. Recently, it has also been used in fractal analysis of bifurcations. 
We study fractal properties of the unit-time map near nilpotent nonmon- 
odromic singularities of planar vector fields using normal forms. We are 
interested in nilpotent singularities because they are nonhyperbolic, and 
we know that near nonhyperbolic singularities the box dimension is non- 
trivial. We study discrete orbits generated by the unit-time map, on the 
separatrices at the bifurcation point, and get results for the box dimen- 
sions of these orbits. Box dimension results will be illustrated in details 
using the examples of Bogdanov-Takens bifurcation and bifurcation of 
the nilpotent saddle. The study is also been extended to the appropriate 
singular points at infinity of normal form for the nilpotent singularity. 
Moreover, we study the unit-time map of the normal form for a saddle, 
near singular points at infinity. 
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1 Introduction 

Since the 1970s fractal analysis has become a tool for studying dynamical sys- 
tems, mostly their invariant sets and related measures. Complexity of invari- 
ant sets, measures, and graphs of functions have been studied using numerous 
fractal dimensions. We stress widely used Hausdorff dimension and box dimen- 
sion (also known as box counting, Minkowski dimension, Minkowski-Bouligand 
dimension, capacity dimension, limit capacity). This approach generated vari- 
ous results concerning box and Hausdorff dimensions of strange attractors, like 
Lorentz or Henon, and also results about Smale horseshoe, Julia and Mandelbrot 
sets. Results about homoclinic bifurcations and fractal dimensions have been 
obtained, also. See survey article |33| and references therein, to find mentioned 
results. 

The new approach of studying dynamical systems by using the fractal dimen- 
sion showed up recently. The idea is quite simple, to compute box dimension of 
any trajectory, and to connect the obtained result with some other properties of 
the studied system. Results about box dimension and Minkowski content can 
be found in recent development of application of fractal analysis to solutions of 
differential equations and dynamical systems. See for example [T3], [H], [T8] . 
[19], [29]. For our study of discrete systems it is particularly interesting frac- 
tal analysis of bifurcations of discrete dynamical systems (see [4], [8], [9]), [27] . 
[3"T] . [TO] , also applied to continuous systems. Direct connection between the 
box dimension of trajectories of dynamical systems and the bifurcation of that 
system has been proved in the articles. 

The first article in that direction was inspired by book of C. Tricot [2"B"] . 
A new insight to fractal dimensions could be found in the book. We learnt 
that nonrectifiability (infinite length) of the curve could be measured by box 
dimension near the point of accumulation. There are two interesting formulas 
in that book, which we exploited very much. Formula for box dimension of 
nonrectifiable spiral r = <p~ a , < a < 1, and formula for box dimension 
of nonrectifiable chirp f{x) = x Q sina; - ^, < a < (3. Tricot's result about 
spiral authors applied to spiral trajectories of planar vector field in article [27] . 
Box dimension of spiral trajectories near weak focus or limit cycle is related to 
the Hopf and Hopf-Takens bifurcation. In that article it is observed that box 
dimension is related to the "potential" of the system to produce limit cycles 
under a small perturbation. This observation made some kind of connection 
between box dimension and famous 16th Hilbert problem, so we think that this 
subject deserves to be explored. The study was extended involving Poincare 
map (first return map) near focus or limit cycle. The Poincare map generates 
one-dimensional discrete system, so results about Poincare map from |31) are 
based on results about discrete systems from [3]. 

Here we are interested in discrete systems, so we stress some results about 
discrete systems. First of all, notice that, for an orbit of one-dimensional discrete 
dynamical system, the Hausdorff dimension fails to show difference between 
systems. Namely, because of its property of countable stability, the Hausdorff 
dimension does not 'see' the countable sets at all. On the other hand, the box 
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dimension is only finitely stable so it can 'see' them clearly. That is the rea- 
son why the box dimension is prefered for studying orbits of discrete dynamical 
systems. It is known that the box dimension of an orbit near hyperbolic fixed 
point of the one-dimensional discrete dynamical systems is trivial (see [J ) . The 
analogous result for the hyperbolic fixed point of the systems in W 1 is showed 
(see [IH])- On the other hand, the box dimension near the nonhyperbolic fixed 
point of discrete dynamical system is strictly positive. The article |8] proved 
the connection between the value of that box dimension and appropriate one 
and two-parameter bifurcations. Also in the article [9] the box dimension result 
for the Neimark-S acker bifurcation is proved. Recent work [35] shows inter- 
esting connection beetween box dimension and Minkowski content with formal 
classification of parabolic diffeomorphisms. 

Discrete systems and continuous systems with spiral trajectories have been 
studied using box dimension. Poincare map makes a connection between con- 
tinuous and discrete systems in the case with spiral trajectories, we call it mon- 
odromic case. For nonmonodromic case the connection could be established 
using unit-time map (also called time 1 map). Here we concentrate on nilpotent 
nonmonodromic singular points, while monodromic nilpotent singularities has 
been studied in [23] . In |21j Hopf-Takens bifurcation at infinity has been studied 
using some generalization of box dimension and Poincare compactification. Ex- 
porting this idea we also investigate unit-time map of nilpotent nonmonodromic 
singularities near singular points at infinity. More about dynamics near nilpo- 
tent singularities and near infinity could be found in pQ, 0, 0, [H], [2S], etc. 

Now we recall the notions of box dimension and Minkowski content. For 
further details see e.g. [5], and for some generalizations see [16], [21]. Let A C 
be bounded. Minkowski sausage of radius e around A is a e-neighborhood 
of A, that is A e = {y G R N : d(y,A) < e}. Let s > 0. The lower and upper 
s- dimensional Minkowski contents of A are defined by 

Mi(A) :=liminf -JM M* a (A) -limsup^L 

Then the lower and upper box dimension are defined by 

dim p A = inf{s > : M%{A) = 0}, dh^i B A = mf{s > : M* S {A) = 0}. 

If dxm B A = dim^A we denote it by dims A. If there exists d > such that < 
■M.i(A) < M.* d (A) < oo, then we say that set A is Minkowski nondegenerate. 
Clearly, then d — dims A. If |^4 e | ~ £ s for e small, then A is Minkowski 
nondegenerate set and dim B A = N - s. If M%{A) = M* S {A) = M d {A) <E 
(0, oo) for some d > 0, then A is said to be Minkowski measurable. Clearly, then 
d = dims A. Let A and B be two disjoint bounded sets such that dims A — 
dims B. It is easy to see, using the finite stability of the upper box dimension 
(dims {A U B) = max{dinisA, dim^S}), the monotonicity of the lower box 
dimension (dhn B (A U B) > dim B A) and dhn B (A U B) < dime (A U B) (for 
details see [5]), that 



dim B (A U B) = dim B A = dim B B. 



(1) 
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In the paper the following definitions are used. We say that any two sequences 
(a„)„>i and (6„)„>i of positive real numbers are comparable and write a n ~ 
b n as n — > oo if A < a n /b n < B for some A, B > and n sufficiently big. 
Analogously, two positive functions /, g : (0,r) — > R are comparable and we 
write f(x) ~ as x — > if f(x)/g(x) € [A, £?] for a; small enough. 
Hence we consider a discrete dynamical system 

x„ +1 = F(x„), Xi e R N 

generated by a C k function F : R N H> M. N . The orbit of a system is a sequence 
(x„)„>i such that x„ + i = F(x„) for some Xi G M. N . Let xo = be a fixed 
point (F(xq) = xo) of that system and let A be a Jacobi matrix DF(xq) at 
xo. The eigenvalues Ai,...,Ajv of the matrix A are called the multiplicators 
of fixed point. We denote by Nq the number of multipliers on the unit circle, 
by ./V_ the number of multipliers inside the unit circle and by N + the number 
of multipliers which lies outside the unit circle. The fixed point is hyperbolic 
if Nq = 0, that is, there is no multipliers on the unit circle. Hyperbolic point 
is called a hyperbolic saddle if N-N + ^ 0. The fixed point is nonhyperbolic if 
Nq^O. 

In this paper the main object of our study is a box dimension of the orbit 
of the unit-time map of planar vector fields near the nilpotent singularity. The 
unit-time map generates a two-dimensional discrete dynamical system which 
correspond to appropriate planar vector field. We know that the unit-time map 
on the characteristic orbit near the hyperbolic singularity is zero, while near the 
nilpotent singularity is positive and we will see that it depend on the order of a 
system and the asymptotic expansion of separatrices. 

The remainder of this paper is organized as follows. In Section 2 we recall 
the possible topological types of nilpotent singularities in the plane, and intro- 
duce the model system which we will study. Then we present the unit-time 
maps for a model system with the nilpotent singularity. In Section 3 we apply 
quasihomogenous blow-up method to the model system, and get the asymp- 
totic behavior of the separatrices, in the cases of nilpotent singularities which 
have the separatrices. In Section 4 we present the main result about the box 
dimension of the orbits of the unit-time map on the separatrices of nilpotent 
singularities. Section 5 gives some examples of the bifurcations in the nilpotent 
singularities such as cusp and saddle. For the cases of nilpotent saddle-node, 
node and the nilpotent singularity with hyperbolic and elliptic sector, we give 
one example of the appropriate system in order to illustrate the results in the 
previous chapters. In Section 6 we study singularities at infinity for nilpotent 
model system, while in Section 7 we give some additional remarks about singu- 
larities which are not nilpotent. With the same technique we study normal form 
for saddle. We introduce a new notion of dual box dimension, which is related 
to dual Lyapunov constants, also called saddle quantities, see [IT] . We show 
an example where the saddle quantities of saddle normal form at the origin are 
related to the unit-time map of singularities of the normal form at infinity. 



Nilpotent singularities 



4 



2 Unit-time map 

First we will recall the known theorem about the classification of nilpotent 
singularities for planar vector fields (see [3]). We omitted trivial cases (1) and 

(2) . 

Theorem 1 [3] (Nilpotent Singular Points) 

Let (0, 0) be an isolated singular point of the vector field X given by 

x = y + A(x,y), 

V = B(x,y), (2) 

where A and B are analytic functions in a neighborhood of the point (0,0), 
and jiA(0,0) = j\B(0,0) — 0. Let y — f{x) be the solution of the equation 
of y + A{x, y) = in a neighborhood of the point (0, 0) and consider F(x) — 
B{x,f{x)) and G(x) — (^ + f{x)). Then the following holds: 

(3) If G(x) = and F{x) = ax m + o(x' n ) for m € N, m > 1, a ^ 0, then 

(i) If m is odd and a > 0, then the origin is a saddle; and ifa<0, then it is a 
center or a focus; 

(ii) If m is even then the origin is a cusp. 

(4) If F(x) = ax m + o(x rn ) and G(x) = bx n + o(x n ), m, n G N, m > 2, n > 1, 
a 7^ 0, b 7^ 0, then we have 

(i) If m is even, and 

(11) m < 2n + 1, then the origin is a cusp. 

(12) m > 2n + 1, then the origin is a saddle-node. 

(ii) If m is odd and a > 0, then the origin is a saddle. 
(Hi) If m is odd, a < and 

( 1111 ) Either m < 2n + 1, or m = 2n + 1 and b 2 + Aa{n + 1) < 0, then the origin 
is a center or a focus. 

(1112) n is odd and either m > 2n + 1, or m — 2n + 1 and b 2 + Aa(n + 1) > ; 
then the phase portrait of the origin consist of one hyperbolic and one elliptic 
sector; 

(1113) n is even and m > 2n + 1 or m = 2n + 1 and b 2 + Aa[n + 1) > 0, then the 
origin is a node, (b > repelling, b < attracting). 

We are interested in nonmonodromic isolated nilpotent singularities, it means 
that the singularities have separatrices or characteristic orbit going through the 
nilpotent singularity. The types are: cusp, saddle, saddle-node and node, and a 
nilpotent singularity with one elliptic and one hyperbolic sector. 

We consider the model system: 

x = y 

V = f(x)+yg(x) + y 2 B(x,y), (3) 

where /, g and B(x,y) are C°° functions, ji/(0) = 5(0) = jooB(0, 0) = 0. Also 
it is satisfied joo/(0) 7^ 0, moreover f(x) — ax m + o(x m ), a^0. For g there are 
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two cases: ioofl'(O) = or g(x) = bx n + o(x n ), b ^ 0. So we consider the system 
x = y 

y = ax m + bx n y + y 2 B(x,y)+o(x m ) + yo{x n ) (4) 

under the assumption deg(_B) + 2 > max{m, n + 1}. 

Before proving the theorems, we recall the procedure for calculating the unit- 
time map of continuous system by using the Picard iterations. So, we consider 
the continuous dynamical system 

x = F(x) (5) 

where x £ M. N , F : l w — > M. N . The simpliest way of getting the discrete 
dynamical system from the continuous one is by using the unit-time map ^t(x). 
Namely, we fix to > and we consider the system which is generated by the 
iteration of the map (f> tQ (map with displacement t along the trajectory of 
(|5|). If we take t = 1, we get the discrete dynamical system generated by the 
unit-time map 

xm-^(x). (6) 

It can be easily shown that the isolated fixed points of ^ corresponds to the 
isolated singularities of ([5|. In order to study the connection between the hy- 
pcrbolicity and stability of these points, we need to find the connection between 
the corresponding eigenvalues of DF(x ) and -D0i(x o ). 

Now we look at the continuous dynamical system with the singularity xo = 

x = F(x) = Ax + F^(x)+F< 3 )(x) + ..., xel", (7) 

where A = DF(0) and F^) are smooth polinomial vector function of order k: 
F«(x) =C(||x|| fe ) i 

Fi ' ) ( x ) = bij li ... j^x^x^ 2 ■ ■ ■ x 3 n . 

We denote the corresponding flow of Q with 0t(x). Now we would like to find 
Taylor expansion of 0t(x) near Xo = by using the process of Picard iterations. 
Namely, let x^(t) = e At x be a solution of linear equation x = Ax with the 
initial value x, and define 

X (fe+D (t ) = e At x + r* e Mt-T) {F (2) ix (k) {T)) + + F (*+i)( X (k)( T ))) dr . 

Jo 

It is easy to show that (fc + l)-iteration does not change the terms of order I < k. 
By the substitution t = 1 in x( fc ) (t) we get the Taylor expansion of the unit-time 
map 4>i ( x ) until the terms of order k 



<Mx) - e A x + g( 2 >(x) + . . . + g( fe )(x) + 0(||xf +1 ), 



(8) 
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where are polynomial vector function of a form as functions F' 1 '. We get 
B = D(j)i(0) = e A , where A = DF(0). It means that xo = is a hyperbolic 
(nonhyperbolic) singularity of ^ if and only if Xo = is a hyperbolic (nonhy- 
perbolic) fixed point of map (|6|. In dimension one, it is obvious because e° = 1. 
In the plane, we have three cases: 

• A has two different real eigenvalues Ai 7^ A2 =>■ B has two different real 
eigenvalues e Al and e A2 

• A has one real eigenvalue A =>■ B has one real eigenvalue e 

• A has two complex conjugated eigenvalues X\ t2 = a ± hi => B has two 
complex conjugated eigenvalues e a (cos b ± i sin 6) 

Using the above procedure, we can prove the following lemma. 
Lemma 1 (The unit-time map) 

Let (xo,yo) — (0,0) be a nilpotent singularity of the system Q. Then the 
following holds: 

1. If m < n + 1 then the unit-time map has a form 

x k+1 = xk + yk + ^xf + ac 11 x™- 1 y k + . . . + ac lm yj? + 0(\\x\\ rn+1 ) 
y k+1 = y k + ax n k l + ad 11 x k n ' 1 y k + . . . + ad lm y^ + 0(\\x\\ m+1 ); (9) 
with the constants cu = cu{m), du = du(m). 

2. If m = n+ 1 then the unit-time map has a form 

x k+1 = x k + y k + ^x k n + ac 21 x m ~ 1 y k + . . . + ac 2 , m ^ 1 x k y k n ~ 1 + 0(||a;|| m+1 ) 
Vk+i = y k + ax k n + ad 21 x k n - 1 y k + ...+ad 2m y k n + 0(\\x\\ m+1 ); (10) 
with the constants c 2 i — c 2 i{m), du = d 2 i(m). 

3. If m > n + 1 then the unit-time map has a form 

x k +i = x k +y k + h -x n k y + bc 3X xl- x y\ + ... + bc^ n _ x x k yl + T^Vk +1 + °(\\ x \\ n+2 ) 

y k+1 = y k + bx n k y k + bd^x^yt + ... + bd 3 ^ lXk y 7 k l + ^V k +1 + 0(\\xf +2 

with the constants csj = csi(n), d%i = d^n), i = 1, . . . ,n — 1. 
Proof. 

In the case m < n+ 1, by using the above procedure we find the Taylor expan- 
sion of the unit-time map up to terms of order m, while in the case m > n + 1 
we can get the Taylor expansion up to terms of order n + 1M 
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3 Separatrices 

We continue with the model system: 
x = y 

V = f(x)+yg{x)+y 2 B(x,y) (12) 

where /, g and B(x,y) are C°° functions, ji/(0) = gt(0) = jooS(0,0) = 0. We 
study the systems with joo/(0) ^ 0, that is f(x) = as™ + o(x m ), a ^ 0. For g 
there are two cases: joog(0) = or g(x) — bx n + o(x n ), b ^ 0. 

In order to find the asymptotics of the separatrices we use quasihomogenous 
blow-up (see [3J, [32 )■ In the quasihomogenous blowing-up procedure we make 
a finite sequence of changes of variables in the system, leading to the desingu- 
larisation of the system. We choose new variables with respect to the Newton 
diagram of the system, it depends on the slope of the side in the diagram. 
For nilpotent singularities it is possible to desingularise vector field using only 
one quasihomogenous blow-up. Regarding Newton diagram, we have 3 different 
cases for nilpotent singularities: 

1. Hamiltonian like case (m < 2n + 1) 

2. Singular like case (m > 2n + 1) 

3. Mixed case m = 2n + 1 

Now we would like to find the asymptotics of the separatrices for each case. 
1. Case m < 2n + 1 



The system ( 12 ) is C°° is equivalent to 



x = y 

y = <5.T m +y(fa;™ + o(x n )) + 0(y 2 ). (13) 

Case 1A m even: 

For 6=1 we have a cusp. Using the homogeneous blow-up: 

x = u 2 



u m+1 y (14) 



we get the system 



1 

~ u y 



+ (15) 



with the singularities: (u, y) = (0, ±y ^fj), which are hyperbolic saddles. Now 
we find the Taylor approximations for the stable and unstable manifolds of this 
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systems, in order to find the approximation for the separatrices of the initial 
system. We look at the unstable manifold for the singularity (0, \J~^pi), and 

the stable one for (0, — yj ^qrj)- The invariant unstable manifold is 



y=\l — — r + ctiu + ot 2 u 2 + 0(u 3 ). 
to -\- 1 



The stable manifold is 



V = -\ — T + Piu + fW + 0{u 6 ), 
m + 1 



and the corresponding separatrices are: 
unstable 



l m + l m + 2 m+3 m+4 

y=\ -x ^ +aix 2 +a 2 x 2 +0{x ? ); 

" m + 1 



stable 



V2 m + l ™ + 2 + 3 ™+4 
—x ^ +/3 lX 2 + /3 2X 2 +0(x 2 ). 
TO + 1 



to + : 

So, the asymptotic behavior of the separatrices near the origin (cusp) is 



2 m + l 

This is what we need in order to calculate the box dimension of the unit time 
map. 

Remark 1. Since the asymptotics depends only on to, it is interesting to 
explore where is hidden the power n. In fact, it is easy to see that the second 
term in the asymptotic scries of scparatrix depend on n. If we include the 
separatrices in the system, we get that a\ = . . . = a^-i = and ^ for 

k = 2n + 1 — m. For example, if m = 2 and n = 1, we get y = \f\x'i + a\x 2 + 

OLixi +. . ., while for m = 2 and n — 3, we have y — \j\xi +a$x i +a w x l t +. . .. 
Case IB to odd: 

For 8 = 1, we have a nilpotent saddle. The case 6 = — 1, where we have a center 
or focus, we do not study here. Using the homogenuous blow-up: 

x = u 

y = u^y (16) 

we get the system 

it = uy 

y - (i-^±f) + o(u), (i?) 
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with the singularities: (u,y) — (0, ±y n 2 +1 ), which are hyperbolic saddles. 

Then, we find Taylor approximations for the stable and unstable manifolds 
of this systems, in order to find the approximation for the separatrices of initial 



system. We look at the unstable manifold for the singularity (0, y ^xj)i an d 



the stable one for (0, —<J ^p±)- The invariant unstable manifold is 

~2 



y = \ j + a\u + a 2 u 2 + 0(u 3 ), 

" to + 1 



and the stable manifold 



'/ = -\l + + P 2 u 2 + 0(u 3 ). 

" m + 1 



Corresponding separatrices are: 
unstable 



I m + l m+3 m + 5 ™ + 7 

J/=y 72: 2 + aix 2 + a 2 x 2 +0(x 2 ); 

-~, + 1 



stable 



y = -\ t x 2 + Pix 2 + fox 2 +o(x 2 ). 

m + l 



So, the asymptotic behaviiour of the separatrices near the nilpotent saddle is 



Z m+l 

w ~ ±W x 2 . 

- + 1 

Notice that it is the same as in the cusp case. 

Remark 2. It is easy to see that the second term in the asymptotic series 
depend on n. If we include the separatrices in the system, we get that the 
separatices are: 



l m+ i m+l j , m+l , ou 

y = ±\ —rx 2 +a k x 2 +L +a 2k x 2 +2A + . . . 

V TO + 1 

where k = 2n +i-' m , p r example, for to = 3 and n = 2, we get i/ = \J\x 2 + 

a>ix 3 + a 2 x 4 + . . .. 

2. Case m > 2n + 1: 

The system ( 12 ) is C°° equivalent to 

x = y 

y = ax m +y(x n +o(x n )) + 0(y 2 ), a^O (18) 

Nilpotent singularities in this case are nilpotent saddle, saddle-node, node and 
singularity with the elliptic and hyerbolic sector. Using the homogenuous blow- 
up: 

x — u 

V = u n+1 y (19) 
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we get the system 

u = uy 

y = y(l - (n + l)y + 0(u)) + an" 1 - 2 "- 1 (20) 

with two singularities: T\(u, y) — (0, ^xx) and T 2 (u, y) = (0, 0). The singularity 
Ti is a hyperbolic saddle, while T 2 is a semi-hyperbolic point with the unstable 
manifold u = 0, and the center manifold transverse to it. The semi-hyperbolic 
singularity T 2 can be saddle, node or saddle-node. Topological type of the 
singularity of the system (18) depends of the type of T 2 . 

Now, we look at the unstable manifold for the singularity Ti(0, ^zj)- The 
invariant unstable manifold is 

y = h ct\u + a 2 u 2 + 0(u 3 ), 

n + 1 

so the corresponding unstable separatrix has a form 

y = -^—x n+1 + a lX n+2 + a 2 x n+3 + 0(x n+i ). 
n + 1 

Then, the asymptotic behavior of the "upper" separatrix near the origin is 

y~ — x n+1 . 
y n + 1 

This is what we need in order to calculate the box dimension of the unit-time 
map. 

Remark 3. Since the asymptotics depends only on n, it is interesting to 
explore where is hidden the power m. In fact, it is easy to see that the second 
term in the asymptotic series depend on m. If we include the separatrices in the 
system, we get that ot\ = . . . = otk-i = and — , for k = 2n + 1 — m, 
that is, 

y— ~~r x + O-m- {2n+\)X + U(X . 

Now, we find the center manifold for the singularity T 2 — (0,0). We know 
that has a form 

V = l\u + 72« 2 + 0(u 3 ), 
so the corresponding invariant manifold for the initial system is 

y = 7ia;" +2 + j 2 x n+3 + 0{x n+A ). 

But, if we include it in the system, we get that the first term also depends on 
m. Since m > 2n + 1, there exists k > such that k = m — 2n — 1. We obtain 
the asymptotic series 

y = i k x n+k+1 + o(x n+k + 2 ) = i m - 2n - lX m - n + o(x m - n+1 ). 
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3. Case m = In + 1: 

The system ( 12 ) is C°° is equivalent to 

x = y 

y = ax 2n+1 +y{x n +o(x n )) + 0(y 2 ), a^O. (21) 

In this case we are interested in the nilpotent singularities such as saddle, node 
and singularity with the elliptic and hyperbolic sector. 
Using the homogcnuous blow-up: 

x = u 

V = u n+1 y (22) 

we get the system 

u = uy 

y- = a + y-(n+l)f + yO(u) (23) 
with several possibilities: 

1. For a > 0, we have two singularities — (0, ^^(n+t)^ 1 ) wmcri are 
both hyperbolic saddles, and the corresponding topological type of the 

initial system is a nilpotent saddle. We denote by A 12 — 1± " V 2(k+i) 1 ^ +1 1 
where A\ > and A 2 < 0. Then the unstable separatrix is of a form 

y = A lX n+1 + a 1 x n+2 + a 2 x n+3 + 0{x n+A ), 

and the stable is 

y = A 2 x n+1 + p lX n+2 + p 2 x n+3 + 0(x n+4 ). 

2. For a < 0, we can have two singularities T\^ 2 = (0, ~^g^^y^^ )i where 
T\ is a hyperbolic saddle and T 2 is a node; or if 4a(n + 1) + 1 = 0), we 
have only one singularity T\ — (0, 2 (n+i) ) wm ch is a saddle- node. By 
blowing down, we can get the nilpotent node (for n even) or a nilpotent 
singularity with one elliptic and one hyperbolic sector (for n odd). The 
invariant unstable separatrix for the singularity T\ is: 

if 1 + 4a(n + 1) > 0, then 

y = A lX n+1 + ai x n+2 + a 2 x n+3 + 0{x n+A )- 
if 1 + 4a(n + 1) = 0, then 

-1 n -L1 

V 



2(n+l)' 



3. If Aa(n + 1) + 1 < 0, then there is no singularities, and the initial system 
has a center or focus. 

We showed that in all nonmonodromic cases y ~ x n+1 . 
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4 Box dimension of the unit-time map 

As we mentioned before, it is already known that the box dimension of each 
orbit of discrete dynamical system near the hyperbolic fixed point in M.-^ is 0. 
Also, it is known that the box dimension of orbit near the nonhyperbolic fixed 
point of discrete system in R w is strictly positive, see [TU]. These results can 
be applied to the unit-time map of the continuous system near nonhyperbolic 
singularity. 

Theorem 2 |10j Let (0, 0) be a hyperbolic singular point of continuous planar 
dynamical system. Then the unit-time map on each characteristic trajectory 
near (0,0) has positive box dimension. 

The results of box dimension in the case of the nonhyperbolic singularity 
with only one multiplier on the unit circle can be found in |10j . Now we would 
like to get the result for the exact value of box dimension near a nilpotent 
singularity. 

Lemma 2 Let A = {xk}ken an d B = {yu}keN be a two decreasing sequences 
which tends to with initial points xq and yg and with the properties x^ — Xk+i — 
x%, for a > 1 and y k - y k+1 ~ y%, for f3 > 1. Let S(x ,y ) = {{x k ,yk)} be 
a two-dimensional discrete dynamical system, with initial point (xo,yo). Then 
the following holds: 

(i) if a > (3. then dinig S = 1 — A; 

(ii) if a < j3 , then dims S = 1 — \ ■ 

Proof. 

From Theorem 1, [3], it follows that dims A = 1 — A and dims B = 1— jj. 
It is obvious that the set A is an orthogonal projection of the set S on the 
x-axis. Analogously, the set B is an orthogonal projection of S on the y-axis. 
Since the orthogonal projection is a Lipscitz map, then we know that dim ^S 1 > 
max{dimB A, dime B}. That is the lower bound for the box dimension of the 
set S. The upper bound for the box dimension can be calculated directly by 
estimating the area of the Minkowski sausage of radius e. We denote by ua{s) the 
minimal n £ N for which x n — x n+ \ < 2s, and analogously for ub(s). Also, we 
denote by ng (s) the minimal n € N such that \]{x n — £ n +i) 2 + (y n — Vn+i) 2 < 
2e. It is easily seen that in the case a > f3, it holds 

n A (V2e) < n s {e) < n B (V2e). (24) 

In the case a < ft, the opposite inequalities are valid. Now, we can estimate 
the area of the Minkowski sausage of radius e which we divide into the tail 
\S e \ t (before overlapping), and nucleus \S e \ (after overlapping). We want to 
calculate the upper bound for the box dimension. So we have 



\S E \ t = 7re 2 (n s (e)-l), 
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and 

2 /"'•'• ' 



\S E \ n <7re 2 + ( 9 (x) + S-( J g{x)-S))dx 

Jo 

where y = g{x) is the curve on which the discrete orbits lies, and it can be easily 
seen that g(x) ~ a; 7 , where 7 = ^Ey- Constant S = max n ^S Xn is a maximum 
5 (see Figure 1) such that almost the whole S e (without two semicircles) is 
between the curves g{x) + 5 and g(x) — S. In the case a > /?, the curve y = g(x) 
is increasing, and the derivative is also increasing so S = S Xns{e) . 
Now we have 

\S e \ < \S £ \ n + \S e \ t = ne 2 n s (e) + 28x ns{e) . 

Notice that 5 Xn > e and that limn^oo S Xrl — e. So, we can choose e small enough 
such that 5 < 2e. Then we have 

\S B \ < Tre 2 n s {e) + 4ex ns(E) . (25) 



In the case a > j3, from (24) it follows 

x n s (e) < X n A (V2e)J 

that is, 

nsieY^ <n A (V2e)'^- 



Now we put x nA (^/2 £ j ~ riA{£) and uaIs) ~ e ^ a' in the inequality (25), 
divide by e 2 ~ s , and get 

^ <^e s C l£ ^^ +2C 2 e s - 1 n A (e)-^ < 
e s 

So it follows that 



dimes' < 1 ■ 

a 



and we proved the lemma. 



Remark 4. This lemma can also be easily proven by using the Lemma 4 
from [30]. 

Remark 5. In general, y = g{x) from the previous proof can be a piecewise 
spline which numerically aproximates the curve. Also, notice that the discrete 
orbits studied here lie on the curves, that is, trajectories or separatrices of the 
continuous system. 
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Figure 1 upper bound for \S e \, a > f3 



Theorem 3 (Box dimension near nilpotent singularity) 

Let we have a system Q). Let (xo,yo) — {0,0) be a nilpotent singularity and let 
y(x) = x 1 + o(x J ), 7 <E (1,to) be a separatrix of the system Q). Let V be an orbit 
on the separatrix generated by the unit-time map of the system near (0, 0), and 
T x , T y are projections ofT to the coordinate axes. Then the following holds: 
(1) If m < n + 1, then dim s = 1 — A dims r y = 1 — ^ . 

(1) Ifj 2 > TO, then dim^ T = dim^ r x = 1 — A. 
(ii) If J 2 < m, then dim# T = dim# T v = 1 — 

(2) Ifm>n+l, then dim B T x = l- ±, dim s T y = 1 - 

(%> #7 > |(1 + \/l + 4n), tfien dim B r = dimsT* = 1 - A 
(M,) Jf 7 < |(1 + \/T+ 4n), t/ien dims T = dims r y = 1 — ■ 

Proof. 

(1) Case m < n + 1: 

From Lemma 1, it follows that the asymptotics of the unit-time map are 

and m 

2/fc ~ 2/fe+l - 2/fe 7 ■ 
Using Theorem 1 from [3], we have 

dimsT^ = 1 - 

7 

and 

dims T y = 1 

We denote by a = 7, /3 = — , then we have 

a > p <^4> 7 2 > m 
a < /3 7 2 < m. 



7^ 
m 
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So the results for the box dimension follow from Lemma 2. 
(2) Case m > n + 1: 

From Lemma 1, it follows that the asymptotics of the unit-time map are 

7 

Xk ~ x k+l — X k , 

and 

- yfe+i - Vk ■ 

Using Theorem 1 from we have 

dim B T x = 1 

7 



and 



dim^ L„ = 1 



n + 7 



Denoting by a = 7, ,3 = + 1, for 7 > we get 

a > /? <^> 7 > -(1 + vT+4n), 
a < f3 <^ 7 < -(1 + \/T+ 4n). 
Now the results for the box dimension easily follow from Lemma 2. 



The results from the theorem will be illustrated by examples in Section 5. 
Remark 6. It is interesting to notice that for cusp and nilpotent saddle, 
we have a characteristic set of values concerning box dimension. For cusp 

_ dim B r B _ _4^_ . _ f 4 8 12 16 20 

c_i dim B r y ~ 2fc + l fc 1 t3' 5' 7 ' 9 ' 11' 

The set D c coincide to the set of values of box dimensions of the spiral trajectory 
near weak focus, see Theorem 9, [27]. Also for the nilpotent saddle, the possible 
values are: 

».-i£iH-<*i^» >-{!•!•?§•£■ 

The set D s coincide to the set of values of box dimensions of the spiral trajectory 
near limit cycle, see Theorem 10, [2"T] . 

Remark 7. In the nilpotent case the connection between the box dimension 
and the multiplicity of fixed point, or the cyclicity of singularity should be fur- 
ther explored. Cyclicity means maximal number of limit cycles which could be 
obtain from the system under small perturbation. Unfoldings of nilpotent sin- 
gularities show their complex structure after application of blowing-up method. 
Blowing-up the singularity we find polycycles "inside". Roughly speaking poly- 
cycles are separatrices which are "closed curves" passing through singularities, 
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the simplest cases are saddle-loop and two saddle-loop. Poincare map near 
singularity, limit cycle or polycycle is a standard tool for studying cyclicity. 
Problem is that Poincare map near monodromic nilpotent singularity or poly- 
cycle is not analytic. These fractal methods could be adjusted to such cases 
using some other scale to obtain an asymptotic expansions, for Chebyshev scale 
see [T5]. 

5 Examples 

In this section we will present several examples of the nilpotent singularities in 
order to illustrate the results for box dimensions of the unit-time map on the 
separatrices. In the case of cusp and nilpotent saddle, we will present the whole 
unfolding for the appropriate bifurcation, while in other cases we will give only 
the nilpotent situation. In the following examples we also use results from [4] 
and [5] dealing with box dimension of the orbits of discrete one-dimensional 
systems at the bifurcation point. The corresponding values for nondegenerate 
saddle-node, and period doubling bifurcations are |, |, while for hyperbolic 
orbit of node and focus, box dimension is equal to 0. For saddles box dimension 
is computed on the stable and unstable manifold. The semihyperbolic cases 
can be reduced to the center manifold, see [TO]. To complete the study about 
box dimension of the whole unfolding we also use result about discrete Hopf 
bifurcation called Neimark-Sacker appearing in 2-dimensional systems, see [S]. 
Discrete spiral orbit at the bifurcation parameter has box dimension equal to | . 

5.1 Cusp 

Let us consider the normal form for the Bogdanov-Takens bifurcation 

x = y 

y = fo + fox + x 2 - xy, (26) 

where £ t are parameters. We can see the bifurcation diagram of Bogdanov- 
Takens bifurcation at Figure 2, see more details in [T2]. We denote by H the 
negative part of fo axis, because it is a curve where Hopf bifurcation occurs. 

For fo = fo — we have a cusp and we use Theorem |3j case (1), (i), with 
m = 2, n = 1, 7 = = §, and get dims S — dim^ S x — 1/3, dims S y = \ 
on the separatrices. See Figure 3a. For other cases we use the results from [4], 
jS], and [0J. At region 1 there are no singularities. Furthermore, by passing 
through the curve T— a saddle and a node appear, so it is a saddle-node bi- 
furcation curve. On the center manifold we have dims S = dims S x = ^ and 
dims S y = (Figure 3c). Somewhere in the region 2 the node becomes a focus, 
and crossing the curve H a limit cycle is born. So box dimension on the Hopf 
bifurcation curve is | (Figure 3e). Passing through the curve P saddle homo- 
clinic bifurcation occurs, that is a saddle-loop is apeared (for recent result about 
P see [20] )■ In region 4 the saddle-loop is broken and there are two singularities, 
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a saddle and a node. If we continue the journey clockwise and finally return to 
region 1, once more a saddle-node bifurcation occurs (curve T + ). 

All such objects are unfolded in the cusp with dims S = 1/3, for /3\ = = 0. 
Notice that the box dimension is nontrivial when some local bifurcation occurs. 
To detect the global bifurcation on P, box dimension near homoclinic loop 
should be computed, see [IB]. All hyperbolic cases inside the regions have trivial 
box dimensions. At the following figures, the trajectories of continuous system 
are drawn in blue colour, while the orbits of discrete system generated by the 
unit-time map are drawn in red. The separatrices are green. Accumulation of 
red points on the separatrix has been measured by the box dimension. 

ft 




Figure 2 bifurcation diagram 




Figure 3c curve T-, dim^ S = 1/2 Figure 3d region 2 from Figure 2 
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Figure 3e curve H, dims S — 4/3 Figure 3f region 3 from Figure 2 




Figure 3g curve P 



Figure 3h region 4 from Figure 2 




Figure 3i curve T+, dims S = 1/2 
5.2 Nilpotent saddle 

The normal form for the two parameter bifurcation of nilpotent saddle is 



x = y 

y = fax + (3 2 y + x 3 - x 2 y, 



(27) 



where /3i.2 € K are parameters. We can see the bifurcation diagram for the 
degenerate Bogdanov-Takens bifurcation at Figure 4, see more details in [12j . 
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For /3i = /?2 = we use Theorem^! case (1), (i), with m = 3, n = 2, 7 = = 
2, we get dims S = dims 5 X = 1/2 and dims S y = 3. See Figure 5a. For other 
cases we use results from [4], [8], and [9]. 

ft 



A 






7- 




H 






(0,0) 




















7- 





ft 



Figure 4 bifurcation diagram 



Now through the bifurcation diagram of the unfolding we can see the chang- 
ing of the box dimension. Similarly as for the cusp case, we start from region 
1 with a hyperbolic saddle (Figure 5b), passing through curve T~ we get two 
more singularities. So in region 2 we have two saddles and one node (Figure 5d). 
This bifurcation is seen by changed box dimension on the curve T~ because on 
one separatrix (orange) is dims S = dims S x — ^ and dims S y — (Figure 5c). 
Then on the curve H Hopf bifurcation occurs (Figure 5e) , while on P two-saddle 
loop appears (Figure 5g). The conclusion is analogous as in the previous case, 
the box dimension is changed at the bifurcation point. Notice that, in this case, 
the box dimension for /3i — fa = is bigger. It could be connected to the fact 
that more object is "hidden" in the nilpotent saddle (3 singularities and 1 limit 
cycle), comparing with the cusp case (2 singularities and 1 limit cycle). 




Figure 5a 0i = 02 = 0, dims S — 1/2 Figure 5b region 1 on Figure 4 
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5.3 Nilpotent saddle-node 

The system 

x = y 

V = x A + xy (28) 

is singular like case m > 2n+l, m = 4, n = 1, from Section 3, and the topological 
type is a nilpotent saddle-node, see Theorem [I] (i2). We use Theorem |3] to 
obtain the box dimension. For this case, we have behavior y ~ x^~ =12 011 
the separatrix, and using Theorem [3] (2), (i) for 7 = 5/2 we obtain dim^ S = 
dims S x = 1 — § = I and dim^ S y = |, where 5 is a discrete orbit generated 
by the unit-time map on the separatrix. See Figure 6. 




Figure 6 nilpotent saddle-node, dims S = 3/5 



5.4 Nilpotent node 

The system 



x 

y 



4x 2 y 



(29) 

is mixed case m = 2n+l, m = 5, n = 2, from Section 3, and the topological 
type is nilpotent node, see Theorem [l] (iii3). We use Theorem [3] to obtain box 
dimension. For this case, we have behavior y ~ x n+1 = x 3 on the separatrix, 



and using Theorem 3 (1), (i) for 7 
and dim^ S y = | 
separatrix. See Figure 7 



3 we obtain dims S = dims S x = 1- 



|, wnere S is discrete orbit generated by unit-time map on the 
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Figure 7 nilpotent node, dims 5 = 2/3 

5.5 Hyperbolic and elliptic sector 

The system 

x = y 

y = -x 3 + 3xy (30) 

is mixed case ro = 2n + 1, m = 3, n = 1, from Section 3, and the topological 
type is a singular point with one hyperbolic and one elliptic sector, see Theorem 
[T| (iii2). We use Theoremkmo obtain the box dimension. For this case, we have 
the behavior y ~ x n+1 — ar on the separatrix, and using Theorem [3] (1), (i) for 
7 = 2 we obtain dims 5 = dims S x = 1 — \ = | and dims 5^ = |, where S is 
a discrete orbit generated by the unit-time map on the separatrix. See Figure 
8. 




Figure 8 hyperbolic and elliptic sector, dims 5=1/2 

6 Singularities at infinity for normal form of nilpo- 
tent singularity 

The Poincare compactification (see [3]) is a standard tool for studying the sin- 
gularities at infinity. In the Poincare compactification, the infinity is presented 
by a circle. 

We continue with the system having a nilpotent singularity at (0,0) 
x = y 

y = ax m + bx n y, a,b^0. (31) 

The case m < n + 1 : 

Using the formulas from [3], p. 152, where 



1 u 

= -, y = -, 

V V 
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and dividing by common divisor, in the first chart we obtain 
u = bu + av n+1 - m - u 2 v n 

y = -uv n+1 . (32) 

Box dimension of discrete system generated by the unit-time map on the center 
manifold v ~ [ s d[m B S(xi,yi) = 1 — 2 n-m+2 > wnere ipXiVx) is a fi rs t 

iteration. 



Figures 9 represent the orbits generated by unit-time map of the system ( 32 1 
with m = n = 2 (Figure 9a) and m = 2, n — 3 (Figure 9b). 





Figure 9b dims S 
Using formulas from [3], p. 152, where 

1 



u 

V 



y = 



and dividing by common divisor, in the second chart we obtain 

u = v n - au m+1 v n+1 - m + bu n+1 

y = 



-au m v n - m+2 



bu n 



(33) 



Box dimension of the unit-time map on the separatrix v 



is dim B 5(xi, yi ) 



1 — ^xj. The result is obtained by computation of the Picard iterations. Figure 
10 represents the unit-time map of the system (33) with m = n = 2 (Figure 
10a) and m — 2, n — 3 (Figure 10b). 
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Figure 10a dims S = § 



Figure 10b dim B S = § 



The case m — n+ 1 has a singularity in the second chart, and the box dimension 
is equal as in the previous case. 

The case m > n + 1 has also a singularity in the second chart, and box 
dimension of the unit-time map on the separatrix v ~ w™- 1 is dims S(xx, y\) = 

m+1 

Remark 8. Let us consider the system (31 1 with cusp at origin (m = n — 2) 



x = y 

y = x 2 + x 2 y, a 7 b^0. (34) 

Our result show that the box dimension near cusp is dimes' = |, while box 
dimensions of cusp near infinity in different charts are dim^ S = | and dim^ S = 
|. It can be interesting to find the connection between these values and the 
cyclicity of cusp. 

7 Singularities at infinity for the normal form of 
a saddle and dual box dimension 

In this section we study the normal form for a saddle, see e.g. [TT], [T5]. In 
[27] and [31], the relation between Lyapunov constants of weak focus and box 
dimension of spiral trajectory has been established, using the normal form and 
the Poincare map, respectively. For a saddle there is analogous normal form 
with coefficients called dual Lyapunov constants or saddle quantities. Here we 
define dual box dimension using connection between weak focus and saddle, 
Lyapunov constants and dual Lyapunov constants. 

Since the saddle is hyperbolic, the unit-time map is exponential and we 
cannot see any interesting behavior of box dimension on the separatrices in the 
singularity. Situation is different at infinity, we show an example where we 
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can see interesting behavior of box dimension depending on the dual Lyapunov 
constants. 

We consider the normal form for a saddle 



x = a x + b y-\ h (x 2 - y 2 ) n (a n x + b n y) 

y = b Q x + a Q y^ h {x 2 - y 2 ) n (b n x + a n y), a n ^ 



(35) 



where the coefficients are called dual Lyapunov constants or saddle quantities. 
If aj. 7^ is fist nonzero coefficient, we say that a saddle is weak of order k, 
completely analogous to the standard definition of a weak focus of order k. 
Saddle quantities play important role in the problem of cyclicity near a saddle 
loop, they appear in asymptotic expansion of the Dulac map. Poincare map 
near a saddle loop is a composition of Dulac map (passing near a saddle) , and 
the regular part. 

Standard normal form for a weak focus is (see |25| ) 



x = —y + a^x + aix(x 2 + y 2 ) ■ ■ • + a n (x 2 + y 2 ) n 

y = x + a y + aiy(x 2 + y 2 ) h a n (x 2 + y 2 ) n , a n ^ 0. 

Since a n ^ 0, we write a n — 1 and in polar coordinates we have 

r = r{r 2n 



(36) 



n-l 
i=0 



(37) 



<p = 1. 



In |27j we exploited Tricot's formula for the box dimension of a spiral T defined 
by r = (f~ a , < a < 1, dim^r = jj^- We obtained the result that if 
cifc 7^ is a first nonzero coefficient then spiral trajectory T has box dimension 
dims r 



4fc 
2fe+l 1 



Put a n = bo = 1, bi = for i = 1 . . . n in (35 1, and change of variables 



r cosh cp, y = r sinh ip 



we get exactly the same system as (37 1, but in hyperbolic coordinates. Solutions 
of system (37) are comparable to r — (p~ a , < a < 1, with appropriate a. If 



the trajectory of system (37) in hyperbolic coordinates is comparable to 
r = ip~ 2k j we say that has dual box dimension dim^j Th = 2 k+i • 

Let study some example of the system (35 ) at oo. Assuming that = bo = 
1, and the other coefficients vanish we get 



x = y + (x 2 -y 2 ) k 
V = x + {x 2 -y 2 ) k . 



(38) 



Using formulas from [3] 



u 
v 
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and dividing by common divisor, in the first chart we obtain the system with 
two singularities (±1,0). We translate it to (1,0), but keep the notation u,v, 
and get 

u = -2uv 2k - 1 +u 2 v 2k - 1 

v = -v 2k -uv 2k - {-l) k u k {2 + u) k . (39) 



The system (39) has a singularity (0,0) and the invariant set is u = 0. Using 
Picard iterations we get that the unit-time map is comparable tow — v 2k . Then 
the box dimension of an orbit of the unit-time map on the v-axes is dim^ S(vi) = 
1 — ^r, where (0,«i) is a first iteration. Box dimension increases by k, which is 
the order of the weak saddle. Results are the same for singularity (—1,0), also 
the second chart does not show any new singularities. 

Figure 11a represents the discrete orbits of the unit-time map of the system 



(|39j) with k = 1 

u = —2uv + u 2 v 
v = -v 2 - uv 2 + u{2 + u). (40) 



Figure lib represents only one orbit of the unit-time map of (40 1 on the i>-axis. 




Figure 11a system (40) 



Figure lib dim B 5(0, vi) = 1/2 



Figures 12a represents the discrete orbits of the unit-time map of the system 
([39]) with k = 2 



•j 3 + u 2 v 3 
v = -v 4 - uv 4 -u 2 (2 + U y 



u = -2uv + u v 

,A „,„.4 „,2/r> , „.\2 



(41) 



At Figure 12b we can see only one orbit of the unit-time map of (41 ) on the 
u-axis. 
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